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Abstract. We present a new statistical analysis that combines helioseismology (sound speed, 
surface helium and convective radius) and solar neutrino observations (the 8 B and 7 Be fluxes) 
to place upper limits to the properties of non standard weakly interacting particles. Our 
analysis includes theoretical and observational errors, accounts for tensions between input 
parameters of solar models and can be easily extended to include other observational con¬ 
straints. We present two applications to test the method: the well studied case of axions and 
axion-like particles and the more novel case of low mass hidden photons. For axions we ob¬ 
tain an upper limit at 3 a for the axion-photon coupling constant of g ai < 4.1 • 10~ 10 GeV _1 . 
For hidden photons we obtain the most restrictive upper limit available accross a wide range 
of masses for the product of the kinetic mixing and mass of < 1.8 • 10~ 12 eV at 3cr. Both 
cases improve the previous solar constraints based on the Standard Solar Models showing 
the power of using a global statistical approach. 
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1 Introduction and motivation 

The Standard Model of Particles (SM) is not a complete theory, and extensions are necessary 
to address some of the most pressing open questions in fundamental physics. Among others, 
some of these are: nature of dark matter, matter-antimatter asymmetry in the Universe, 
origin of neutrino masses, strong CP-violation problem. In order to solve these problems, 
physics beyond the SM is needed and the existence of new particles and/or non-standard 
properties of known particles are generally invoked. 

Stellar physics, while far from being a closed subject, is a mature discipline and it 
provides an accurate understanding of the internal structure of stars and their evolution. 
The sheer mass and size, the extreme conditions reigning in stellar interiors and, in many 
cases, the extremely long lifetimes that allow to integrate small effects over a very long time 
make stars appealing as laboratories for particle physics under conditions not reproducible 
anywhere else. 

In particular, the existence of weakly interacting particles beyond the standard model or 
of standard particles with non-standard properties can modify the internal structure and its 
evolution in different ways. For example, weakly interacting massive particles (WIMPs), too 
massive to be created inside stars, can be accreted from the dark matter halo of the Milky Way 
and contribute to the transport of energy inside stars [1-3] and, in the case of self-annihilating 
particles act as a localized energy source. Weakly interacting light particles, on the other 
hand, can be thermally produced in the stellar interiors and can easily escape due to their 
large mean free path and act as energy sinks. The increased rate of energy-loss then changes 
the internal structure of stars and also modifies evolutionary timescales. The production 
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rate of each particle depends differently on stellar conditions. Consequently, the impact of 
different particles on stellar structure will depend on the class of stars and on the evolutionary 
state. Several works have given constraints on the properties of the weakly interacting light 
particles using different stellar objects. Some examples are: globular clusters, where the 
luminosity of stars at the tip of the red giant branch [4-6] and the lifetime of horizontal 
branch stars [7-10] can be used to give constraints on axion (e.g. strength of coupling to 
photons and electrons) and neutrino (magnetic dipole moment) properties; intermediate mass 
stars for which sufficient cooling during the helium-core burning phase can make the so-called 
blue loop disappear, thus contradicting the observation of Cepheids [11]; white dwarf stars, 
where the anomalous cooling induced by axions and axion-like particles shortens their cooling 
timescale affecting both their luminosity function (LF) [12-15] and pulsation properties [16- 
18]; neutron stars, where additional axion energy losses [19, 20] can even lead to observable 
deviations of surface cooling in real time [21]; SN1987A, where the extra cooling of the proto¬ 
neutron core would shorten the neutrino pulse [22-25] and the axions produced could produce 
a prompt gamma ray signal [23, 26]. These ideas have been widely discussed in the literature, 
the most comprehensive reference being the seminal book by G. Raffelt [27], which can be 
complemented with recent topical reviews on QCD axions [28, 29], axion-like particles [30] 
and hidden photons [9]. 

Many studies have focused on using the Sun for setting limits on the properties of 
different types of particles; we review some of them below. The Sun is by far the best- 
known star. The solar structure, revealed by helioseismology and solar neutrinos, is well 
determined, and accurate solar models give us information about the past, present and the 
future of the Sun [31]. While in some cases (e.g. axions) the most restrictive limits are 
not inferred from solar studies, the Sun remains the most useful benchmark for testing and 
validating both stellar models and, as it is partly the case of the present paper, different 
statistical approaches to constrain particle properties. Also, it is important to keep in mind 
that CAST [32] and the forthcoming IAXO [33, 34] are experiments specifically designed to 
detect exotic particles directly from the Sun, so having predictions of expected solar fluxes 
for exotic particles remains an important aspect to consider. 

Solar constraints on particle properties have been generally derived from applying lim¬ 
its to variations of either neutrino fluxes [35, 36] or the sound speed profile derived from 
helioseismology [35]. However, a systematic approach aimed at combining different sources 
of data accounting in detail for observation and theoretical errors is badly missing in the 
literature. Here, we try to supply such a tool. 

The initial goal of this work is to extend the general statistical approach presented 
in [37] to constrain properties of particles (e.g. mass, coupling constant) making the best 
possible use of the available information on the Sun, both observational and theoretical. To 
do so we use the helioseismic data combined with the neutrino fluxes in a statistical approach 
that includes the theoretical and observational uncertainties and takes into account possible 
tensions among data and solar model input parameters. We then derive solar limits for the 
well-studied hadronic axions -to gauge the performance of our statistical approach- and for 
the more novel case of hidden photons for which the Sun sets the most restrictive limits on 
the kinetic mixing parameter for small hidden photon masses, m <eV. 

Axions are light pseudoscalar particles that were introduced by the Peccei-Quinn [38] 
solution to the strong CP problem. They are very light particles and interact with ordi¬ 
nary particles much like neutral mesons ( 7 r°, 77,...) but with coupling strengths vastly weaker, 
and some model dependencies. The most relevant for astrophysics is whether axions cou- 
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pie with electrons with similar strength than to nucleons, because this coupling tends to be 
very efficient to produce axions in stellar environments. The DFSZ model [39] as any axion 
model embedded in a GUT are examples of axion models with tree-level axion coupling. 
Here, the main axion production mechanisms are the ABC proceses: axio-recombination, 
bremsstrahlung and Compton [40]. In the KSVZ model (hadronic axions) [41, 42] this cou¬ 
pling is absent at tree level and the relevant axion coupling is to two photons. Here, the 
Primakoff effect is the mechanism that converts photons into axions in the presence of elec¬ 
tric or magnetic fields. The inverse Primakoff effect is used in helioscopes like CAST and 
SUMICO (and the proposed IAXO) to convert solar axions into detectable X-rays and con¬ 
strain the solar axion flux. 

The phenomenology of axions is extended to axion-like particles (ALPs) in a completely 
straight-forward way. ALPs are also bosons with a two-photon coupling, but their mass and 
interaction strength are in principle unrelated, unlike the axion case. These particles arise 
copiously in string compactifications (a string axiverse [43-45]) , where 0(100) ALPs can be 
generically invoked, although the number of light enough particles to be relevant for stellar 
physics is completely model-dependent. ALPs can be behind some very interesting hints 
recently pointed out in astrophysics like the transparency of the universe to gamma-rays [46] , 
the 3.55 keV line [47-49]and the soft excess of the coma cluster [50]. 

In this paper, we are primarily interested in the case of hadronic axions because the 
Sun is one of the most sensitive environments to look for their effects. Solar constraints are 
even more important for the case of ALPs, where the SN1987A constraint is absent in generic 
models (ALPs do not generically have large couplings to nucleons but axions do). We shall 
not endeavour constraining the axion-electron coupling because limits from white dwarves 
and red-giant stars are comparatively much stronger. 

Until now, several works have provided constraints to the axion-photon coupling con¬ 
stant g a j using the variations that axions produce on helioseismologic quantities or neutrino 
fluxes. An upper limit for g ai = 10 • 10~ 10 GeV _1 is found by setting a limit to the deviation 
that axions can impart to the solar model sound speed at a given depth in the Sun [35]. This 
work also gives values for the solar neutrino fluxes depending on axion emission and [36] uses 
this relation to give a more restrictive constraint of g (irf = 7 • 10 _10 GeV _1 at a 3<r confidence 
level using the observed 4>( S B) solar neutrino flux measured by the SNO experiment (see [51] 
for the global analysis of the three SNO phases). In [52] they construct the so-called seismic 
models (non-evolutionary solar models constructed in ad-hoc manner to reproduce the sound 
speed derived from helioseismology), with different values for the axion-coupling constant 
and obtain an upper limit of g ai = 2.5 • 10 -10 GeV -1 by comparing the predicted <h( 8 B) 
with the experimental result and using la uncertainties. For the mass range m a ^ 0.02 eV, 
the most restrictive limit comes from the helioscope CAST with g ai < 0.88 • 10 _10 GeV^ 1 
[53]. The future helioscope IAXO should improve these results, as it is expected to reach 
sensitivities to the axion-coupling constant 1 to 1.5 orders of magnitude better than CAST 

[33]- 

Hidden photons, our second case of study, are vector bosons that couple weakly via 
kinetic mixing with standard photons. The kinetic mixing is represented by the parameter x 
and together with the hidden photon mass, m , are the parameters that must be constrained 
(see [9, 54-56] and references therein). In fact, as described in Section 2.3, solar constraints 
are mostly sensitive to the product x m if rn < eV. 

Hidden photons can only be produced from photon -H- HP oscillations, which are affected 
by the photon refraction in the solar plasma. The oscillations are resonant when the hidden 
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photon and photon dispersion relation match and this happens differently for transversely 
polarized photons and longitudinal excitations (L-plasmons). Resonant emission in the Sun 
is possible for hidden photon masses below ~0.3 keV, the highest plasma frequency u> p , which 
is then also the highest produced photon effective mass. The emission of L-modes is more 
important for low hidden photon masses (below m ~eV) [8, 9], for which resonant conditions 
happen all through the solar interior (each region of the Sun emits L-HPs with energy equal 
to the local plasma frequency). This is the case of interest in this paper. Resonant emission 
of T-modes dominates the energy loss in the hidden photon mass range ~ eV-0.3 keV and it 
is localized in a narrow spherical shell of the solar interior for which the hidden photon mass 
matches the plasma frequency m ~ <u p . This case is very interesting too, but our solar model 
codes need to be tuned for such highly localized sources, a task that we plan to endeavor 
in a forthcoming publication. For even higher masses, arguments from horizontal branch 
stars can give better constraints because, owing to the higher plasma frequencies, resonant 
production extends to higher masses, and when it is not possible the higher temperatures 
ensure the existence of photons of sufficient energy to account at least for the hidden photon 
mass. Therefore, it is not crucial to reexamine this regime. 

This paper is structured in the following way. In Section 2 we present the theoretical 
and experimental aspects of this work: solar models and data used and their respective 
treatment of uncertainties. In Section 3 we introduce the statistical method, in Section 4 
we describe our main results and present new upper bounds to the axion mass and kinetic 
coupling of hidden photons. In Section 5 we compare our work to previous results and draw 
some conclusions and a short summary in Section 6. 


2 Models and Data 
2.1 Standard Solar Models 

In this work we use standard solar models (SSMs) as reference models. SSMs are calibrated to 
match the present-day solar radius R 0 = 6.9598-10 10 cm, luminosity L 0 = 3.8418-10 33 erg s~ 3 
and surface metal-to-hydrogen ratio (Z/X) Q . The choice of this last constraint is critical 
because it essentially determines the distribution of metals in the entire solar structure and 
it has been the subject of much discussion over recent years in the context of the solar 
abundance problem [57-60]. 

One frequently adopted solar photospheric composition is that provided by Grevesse 
and collaborators (hereafter GS98) [61], based on 1-D solar atmosphere models, from which 
(Z/X) 0 = 0.0229. SSMs based on GS98 composition are in good agreement with helioseismic 
observables and solar neutrino fluxes determinations. A more recent compilation of the solar 
photospheric abundances has been provided by Asplund and collaborators [62] (hereafter 
AGSS09). It is based on a more accurate analysis with new 3-D hydrodynamical models 
of the solar atmosphere that reproduces very well the observed solar atmosphere properties 
(granulation, convective velocities and limb darkening among others). As a result of these 
and other improvements in solar spectroscopic analysis, they obtain a much lower surface 
metallicity ( Z/X)@ = 0.0178 that leads to SSMs in strong disagreement with the solar 
structure inferred from helioseismology. More quantitatively, SSMs computed with AGSS09 
surface composition disagree with the observed value of the convective radius R 0 z and surface 
helium abundance Yg at 2.5 a and 3.4cr while when using GS98 composition one obtains 
agreement within 0.2<r and 1.1 a, respectively. For the sound speed profile the situation is 
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similar with average rms deviations from helioseismic determinations a factor 4 worse for the 
AGSS09 models compared with the GS98 ones [31, 37]. 

In [63], and more recently in [37], it was pointed out that current seismic and neutrino 
observables are only sensitive to the actual opacity profile 1 that defines the temperature 
stratification of the Sun. The opacity profile is determined by a set of radiative opacity 
calculations (eg. OP [65], OPAL [ 66 ]) and the solar composition (eg. GS98 or AGSS09). 
While it is not possible to separate these two contributions, the solar opacity profile is well 
constrained by helioseismic and solar neutrino experiments [37]. The production of the exotic 
particles considered in this paper and their effects on solar structure depend on the thermal 
structure of the Sun, not on its detailed composition. This is evident by considering that 
energy-loss rates depend mainly on temperature and density. Therefore, the fact that the 
present observational constraints determine very well the thermal stratification of the Sun 
allow us to use these models to test the impact of exotic particles regardless of the detailed 
knowledge of the composition of the Sun. 

We have computed all our models by using GARSTEC (GARching STEllar Code) [67]. 
The physics used in the calculation of SSMs is as follows: the 2005 update to the OPAL 
equation of state [ 68 ], OP opacities [65] complemented at low temperatures with molecular 
opacities [69] and nuclear rates from Solar Fusion II [70]. All models include microscopic 
diffusion as described in [71]. For more details about the SSMs used on this work see [31]. 


2.2 Models with axions 


The production of axions in the Sun occurs via the Primakoff effect, i.e. axions are produced 
by the conversion of photons in the electric field of nuclei and electrons with the interaction 
Lagrangian £ a7 = g ai B • E a, where a is the axion field. Thus, constraints can be placed on 
axion-photon coupling constant g aj . The energy-loss rate per unit mass e a7 is given by [35]: 


3a 7 


c a7 


rnl 

47T p 


( 2 . 1 ) 


where T is the temperature, p the density and F(k 2 ) is a dimensionless function describing 
screening effects given by: 


F(k 2 ) 




( x 2 + K 2 )ln 





Here, the parameter k is defined as: 


( 2 . 2 ) 
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n = ira —3 


Y e + TZ 2 Y 


3 


3 x 3 


(2.3) 


where a is the fine structure constant, ub the baryon density, Y e the electrons per baryon, 
Y j = Xj/Aj and Xj, Aj and Zj represent, respectively, the mass fraction, atomic weight 
and atomic number of the nuclear species j. For solar conditions the function F(n 2 ) can be 
approximated by [35], 

F(k 2 ) = 1.842(k 2 /12) 0 ' 31 . (2.4) 

x The two exceptions are the adiabatic index Id, see [64], that can be used together with the equation of 
state, to determine the metallicity of the solar envelope and the CNO neutrino fluxes that can be used to 
determine the sum of the carbon and nitrogen content of the solar core. 
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Solar models have been computed for different g ar/ = g±o ■ 10 ~ 10 GeV -1 values, where 
gio spans the range from 0 to 20 with an interval of Agio = 1- Two sets of models have been 
computed for the GS98 and AGSS09 solar composition as reference compositions. However, 
following [37], we treat the solar composition as free parameters in our analysis, as discussed 
later in Sect. 3 in more detail. 


2.3 Models with hidden photons 

The production of hidden photons ('y ') in the interior of the Sun can be seen as 7 — 7 ' 
oscillations, described by the Lagrangian 

c = -\a^a^ - + ^-B^ - 1 vr, 

where A^ v and B^ are the held strengths of the photon and hidden photon held, A^ and B^ 
respectively. The energy loss rate per unit mass is calculated by using the approximation 
presented in [9]. We only take into account the dominant process of resonant emission of 
longitudinal hidden photons (details in [9, 54]) for which the energy loss rate is given by 


X 2 m 2 uip 1 
hp e up/T _ 1 4 - 71 - p 


(2.5) 


where x is the kinetic mixing parameter, m the mass of the hidden photon and up is the 
characteristic plasma frequency. Typical values in the solar center are up ~ 0.3 keV and 
T ~ 1 keV. By expanding the exponential, it can be seen that the temperature dependence 
of the energy-loss rate is linear to hrst order in up/T ~ 0.3. In principle, a threshold factor 
0 ^ (m/up) 2 should be included in Eq. (2.5), but this term is completely irrelevant for the 
range of masses considered in this paper, i.e. m < eV. As a consequence, the energy loss rate 
€h p depends only on the product x rn which is the quantity that can be constrained by solar 
data. For hidden photons we consider x m ranging from 0 to 8 • 1CP 12 eV with an interval 
= 1 • ICR 12 eV. As for axions, models with both GS98 and AGSS09 composition as 
reference have been computed, but the complete treatment of composition is given in Sect. 3. 


2.4 Observables 

The observable quantities used in this analysis are the boron and beryllium neutrino fluxes 
4>( 8 B) and 4>( 7 Be), the convective envelope properties inferred for helioseismology, i.e. the 
surface helium abundance Ys and the convective radius RcZ) and the solar sound speed 
profile c a (r). 

The solar neutrino fluxes determined from experimental data are taken from [31], and 
have been derived using all available experimental neutrino data using the method described 
in [72], 

The observational values and errors together with the corresponding SSM predictions 
of the solar neutrino fluxes and convective envelope properties are specified in Table 1. 

For the sound speed profile, we take the 30 radial points 77 in the range 0.05 < t^/Rq < 
0.80 given by [74] and use the helioseismic determination of the sound speed c i in those points 
as our observable quantities. The inversion of the sound speed profile has a dependence on 
the reference model [75] that, however mild, is better to take into consideration. For this 
reason, the solar sound speed profile has been derived using helioseismic data independently 
for each model considered in this work, i.e. using consistently each model for the helioseismic 
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Q 

AGSS09 

GS98 

Observables 

Ref. 

Y S 

Rcz/R© 

$( 7 Be) 

$( 8 B) 

0.232(1 ±0.013) 
0.7238(1 ±0.0033) 
4.56(1 ±0.06) 
4.60(1 ±0.11) 

0.243(1 ±0.013) 
0.7127(1 ±0.0033) 
5.00(1 ±0.06) 
5.58(1 ±0.11) 

0.2485 ± 0.0035 

0.713 ±0.001 
4.82(11™!) 
5.00(1 ±0.03) 

[57] 

[73] 

[31, 72] 
[31, 72] 


Table 1 : The first two columns show the SSM predictions for GS98 and AGSS09 composition and the 
corresponding theoretical uncertainties [37]. Note that theoretical errors do not include uncertainties 
due to the solar composition. The third column summarizes the observational values and errors. 
Neutrino fluxes are in 10 9 cm -2 s -1 for the <f>( 7 Be) and 10 6 cm _2 s _1 for 4>( S B). 


inversion. Inversion has been done using the SOLA inversion technique. The adopted fre¬ 
quencies are the BISON-13 dataset complemented with data of MDI, GOLF and IRIS. More 
details on both the frequency dataset and inversion technique are given in [74]. 


2.5 Uncertainties 


There are two different error sources in this work: the fully correlated theoretical errors, 
induced by uncertainties in the input parameters for solar model construction, and the ob¬ 
servational errors associated to the helioseismic data and the experimental solar neutrino 
fluxes. 

The errors on theoretical predictions are calculated by propagating the uncertainties of 
SSM input parameters I. As in [37], the input parameters are: the age of the Sun, diffusion 
coefficients, luminosity, opacity and astrophysical factors of relevant nuclear reactions (Sn, 
S 33 , S 34 , S 17 , S e 7 , 8144 ). We do not include the chemical composition uncertainty because 
we consider the abundances of volatile and refractory elements as free parameters in the fit, 
as it is explained in Sect. 3. For each observable quantity Q, we calculate the fractional 
change Cqj that is obtained when the fractional variation 51 is applied to the input I. This 
coefficient is obtained as 

Cqj = Bqj 51 (2.6) 


where 


Bqj = 


<91nQ 
(Ill!I 


(2.7) 


The values of the power-law exponents Bqj and of the input parameters uncertainties 51 are 
given in [76] and [37]. The total theoretical error for Q is found by combining in quadrature 
all the input error contributions (see [37] and [59] for details): 


a 


2 _ 

Q,th 


^2 c qj 


( 2 . 8 ) 


Errors for the theoretical neutrino fluxes and convective envelope properties are given in the 
first two columns of Table 1. 

The observational uncertainties Uq are treated as uncorrelated. For neutrino fluxes and 
convective envelope properties, they correspond to the observational errors summarized in 
the third column of Table 1. For the sound speed, a relevant contribution is provided by 
errors in the inversion of helioseismic data. Therefore, at each radial point r t we calculate 
the total uncertainty as: 

Uc(n) = JulpM + ULiri), (2-9) 
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where the first term is the experimental error coming from helioseismic frequencies, while 
the second and dominant term is the so-called ’’statistical” error in the inversion procedure 
estimated by [77]. As in all previous works, we assume that there are no correlations between 
sound speed determinations derived at different radial points. 

Both model and observational sound speed errors are presented in more detail in 
Sect. 4.2. 


3 Method and statistical procedure 

The statistical approach is based upon constructing a y 2 function that can be used a figure- 
of-merit for the quality of different solar models in reproducing the observables. We build 
this function by considering 34 different observable quantities: the neutrino fluxes 4>( 8 B) and 
<f>('Be); the convective envelope properties Y§ and Rcz and the sound speed determinations 
c i = c(ri ) for 30 different value of r/R 0 where r/R 0 < 0.80. We exclude the sound speed 
profile for r/R 0 > 0.80 because the transport is convective in that region, the temperature 
gradient is adiabatic and then, the sound speed profiles of the models are not affected by 
changes in the solar interior and are in agreement with helioseismic results. For values 
close to r/R 0 « 1.0, the predicted sound speed deviates from the observations due to the 
so-called surface or non-adiabatic effects that are not properly accounted for in the models. 
However, solar structure in this region is determined by the surface properties (mainly surface 
temperature and gravity) and, again, not influenced by the properties of exotic particles. 

The fractional difference 5Q between theoretical predictions and observational data is 
defined as 

8Q = 1 - (3.1) 

^cobs 

where Q t h and Q 0 b s denote the predicted and observed values respectively for the quantity 
Q. The differences 8Q are affected by uncorrelated experimental errors Uq and by correlated 
theoretical uncertainties Cqj. Following the formulation of [78], y 2 can be calculated as: 


y 2 = min 
{?/} 


f ~ ^2 1 

AV a Q 


+Y.& 

I 


(3.2) 


where the shifts —£ iCqj describe the effects of correlated errors. They give the corrections of 
theoretical predictions Qth when the input parameters / are varied by the fractional amount 
£i8I. To normalize the effect of these corrections, a penalty Yli £? is introduced in the 
y 2 . The values £/ that minimize the y 2 are referred to as pulls of correlated error sources. 
The quadratic sum y 2 yst = £ 2 give the systematic error contribution to the y 2 while the 

distribution of the pulls gives information about tensions between input parameters [78] in 
the models. We remark that this method is completely equivalent to the use of the covariance 
matrix in order to compute the y 2 function. The main advantage of this approach is that 
it makes it possible to study the individual contributions to y 2 . For more details about the 
statistical method applied to solar models we refer the reader to [37]. 

If we keep the composition of the Sun fixed, there is only one free parameter in our calcu¬ 
lations (either <710 or the product ym) and the energy loss is proportional to its squared value. 
However, we want to prevent that incorrect assumptions about the surface composition of the 
Sun bias our final results. For this reason, we take a conservative approach for constraining 
properties of exotic particles and extend the idea presented in [37] where solar models with 
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arbitrary composition are confronted with observational data. Following [37], we indicate 
with {zj} the heavy element admixture, expressed in terms of the quantities Zj = Zyg/Xg 
where Zj t g is the surface abundance of the j-element, Xg is that of hydrogen, and the index j 
runs over metals. We group the relevant elements for solar model calculations in volatiles (i.e. 
C, N, O and Ne) and refractories (i.e. Mg, Si, S and Fe) and we assume that the abundances 
Zj within each group vary by the same multiplicative factors (1 + <5z vo i) and (1 + fc re f) with 
respect to the values prescribed by the adopted reference solar composition (either GS98 or 
AGSS09). We compute the effects of a composition change by linear expansion in dz vo \ and 
Sz re f with the logarithmic derivatives given in [37] and we superimpose composition effects to 
those produced by axions or hidden photons energy losses. This approach has been shown to 
reproduce well the structural (temperature, density) changes in solar interiors when the solar 
composition is varied [37]. Because axion and hidden photons energy losses do not depend 
on the detailed metal composition of the Sun and the allowed range of energy losses is quite 
limited, this procedure is also valid in solar models including axions or hidden photons. 

In this approach, there are three free parameters in theoretical predictions: one (gio or 
Xm) is related to axions or hidden photon properties and the other two (<5z vo i, Sz Te {) are use d 
to rescale from the reference solar models. Because we are not interested in constraining 
composition, we minimize with respect to (fc vo i, fc re f), i.e. for each value of g\o or %m, we 
choose the solar surface composition that leads to the best agreement with observational 
data. As a result, we obtain the function 

X 2 (/) = min [x 2 if, Sz voh fe ref )] (3.3) 

{fcj} 

where / = g±o, %m, that can be used to constrain axion and hidden photon properties. The 
best fit values for gio and x m are found by minimizing this function and the obtained value 
Xmin provides information on the goodness of the fit. The allowed regions for g\Q and x m 
at 1, 2 and 3er confidence level (CL) are determined by cutting at the values of the variable 
Ay 2 = x 2 — y 2 nin =1,4, and 9 respectively, as it is prescribed for a chi-squared function with 
one degree of freedom. 

As we shall see, the redundancy of the different experimental information included in 
our global analysis allows us to disentangle the peculiar effects introduced by exotic particles 
from those produced by arbitrary composition changes. By using this extremely conservative 
approach, we obtain restrictive bounds which are clearly more robust than those obtained in 
analyses that consider one observable quantity at a time (e.g. the 8 B neutrino flux or one 
point of the sound speed profile, see e.g. [35, 36, 52]). 

4 Results 

Equations 2.1 and 2.5 give the dependence of the energy loss rates induced by axions and 
hidden photons on temperature and density. Different dependences translate into different 
changes in the structure of the Sun. To facilitate discussion and interpretation of results we 
show in Figure 1 the normalized production profile, i.e. the energy loss rate per unit mass 
multiplied by the mass contained in a spherical shell of radius r, as a function of the solar 
radius for axions (red) and hidden photons (blue). The differences can be easily seen. The 
production of axions induces energy-losses in the inner region of the Sun (r < 0.4 Rq), whereas 
the energy-loss distribution is much broader for the case of hidden photons. The difference is 
mostly due to the temperature dependence of the energy-loss rates, see eqs. 2.1 and 2.5. For 
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Figure 1: Normalized production distribution as a function of t/Rq. Red and blue lines correspond 
to axions hidden photons respectively. 


axions, the rate e a7 is roughly proportional to T 6 , while for hidden photons the dependence 
of the rate €h p on the temperature is approximately linear under solar conditions. As a rule 
of thumb, we expect that observables will be affected in a different way. For example, axions 
will have a major effect on the sound speed profile in the core regions while hidden photons 
can produce modifications of the sound speed profile in all regions. Notice, however, that 
although energy losses are localized in certain regions of the Sun, structural changes can be, 
to some degree, present all over the Sun. Indeed, local changes in certain quantities produce 
variations in the whole structure [79]. 

Sections 4.1 and 4.2 show the changes in the solar observables (Sect. 2.4) induced by 
axions and hidden photons in models for which the solar composition is kept fixed at the 
reference values (GS98 or AGSS09). This makes it easy to compare the different response 
of solar models to different particles and to compare with previous work where the solar 
composition was not varied. However, as described before, the solar composition is free in 
the models from which we derive limits on particle properties. The best fit SSM obtained with 
this procedure is briefly presented in Sect. 4.3. Results of this exercise and the constraints 
we derive for both axions and hidden photons are described in 3. 

4.1 Solar neutrinos and convective envelope properties 

Figures 2 and 3 show the dependence of the solar neutrino fluxes 4>( 8 B) and <3?( 7 Be) and 
of the convective envelope properties Yg and Rcz on gio and x m i respectively. Red lines 
correspond to solar models implementing AGSS09 composition and blue ones to GS98. The 
shaded zones depict the la theoretical errors calculated as described in sec. 2.5. The black 
lines shows the experimental values with lcr errors, as given in Table 1. 

In the upper panels of Figure 2, we see how the observables Yg and Rcz change in 
the presence of axions. The surface helium abundance Yg decreases with increasing values 
of the axion-coupling constant because it is almost perfectly correlated to the initial helium 
abundance Yi n i. This quantity decreases with g±o because a higher initial amount of hydrogen 
is necessary to match the solar luminosity Lq in the presence of axion energy losses. The 
change in the convective radius Rcz is quite small because energy losses are localized in the 
innermost regions of the Sun. 
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Figure 2: Evolution of the model parameters (Ys, Rcz, 4>( 7 Be), $( 8 B)) as a function of the axion- 
photon coupling constant. Red color correspond to AGSS09 abundances and the blue one to the GS98 
abundances. Black lines represent the observational value with the errors and the shaded zones show 
model errors. 


The lower panels present results for the neutrino fluxes. <f>( 8 B) and <&(' Be) increase with 
increasing g±o predominantly as a result of higher core temperatures and reach values well 
outside of the theoretical and experimental 1 -er errors already at relatively small values of 
g io- The <3?( 8 B) relative changes are larger than those of <3?( 7 Be) due to its stronger sensitivity 
to temperature. As a consequence, we expect that 8 B neutrino measurements give stronger 
constraints on g\$. It is important to mention that the present experimental determinations of 
d>( 8 B) and <h( 7 Be) have very small uncertainties, 3% and 4.5% respectively. The constraining 
power of both fluxes, and particularly 4 > ( 8 B), is currently limited by uncertainties in solar 
models. 

Figure 2 also illustrates that the relative variations in the observables do not depend on 
the reference solar composition used, as expected because axion production does not depend 
directly on it. Differences between models with GS98 and AGSS09 compositions simply 
reflect the initial differences in the respective SSMs. 

In the case of hidden photons, the changes of the solar neutrino fluxes and convective 
envelope properties have the same qualitative behavior, as it is seen in Figure 3. The decrease 
of Yg and the increase of 4>( 7 Be) and <3?( 8 B) are again a consequence of the solar luminos¬ 
ity constraint, i.e. of the fact that increased initial hydrogen abundance and increased core 
temperature are necessary to compensate energy losses. The behaviour of Rcz is instructive 
because it shows that different types of particles can induce different changes in the solar 
structure. Let us consider, e.g., models for axions and hidden photons for which the boron 
neutrino flux is comparable, say 4>( 8 B) ~ 15 • 10 6 cm 2 s -1 . In the case of axions, this corre¬ 
sponds to solar models with an increase of the convective radius equal to ARcz ~ O.OO3R0. 
For hidden photons, the increase of Rcz is a factor of two larger, i.e. ARcz ~ O.OOGRq. 
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Figure 3: Same as Figure 2 but for the hidden photons case. The variable on the x-axis corresponds 
to the product of the kinetic mixing and the hidden photon mass. 


For a given change in the central conditions, as essentially determined by <f>( 8 B), hidden 
photons lead to larger changes than axions in the outer layers as a consequence of the milder 
temperature dependence of their emission rate. 

A direct comparison between results of Figures 2 and 3 should not be done because they 
are shown as functions of parameters not related to each other in a direct way. To compare 
the effects of both particles, Figure 4 illustrates the fractional change of the neutrino fluxes 
with respect to the SSMs for axions (red lines) and hidden photons (blue lines) as a function 
of the luminosity Lj/Lg carried away either by axions or hidden photons, i.e. Li = L a or 

Lhp- 

For the same Lj/Lg value, changes in the neutrino fluxes are larger for the hidden 
photons models. It was shown by [36] that L a and <h( 8 B) can be related by a simple analytic 
relation of the form 


4>( 8 B) ( L a + L e \ a 

$ssm( 8 B) V L 0 ) 


(4.1) 


where a = 4.6, based on the older generation of SSMs computed by [35]. Our calculations 
yield a = 4.4 for axions, very close to the previous result. Interestingly, for hidden photons 
we find that the same functional form can be used but with a much steeper relation given by 
a = 5.7. 

These results reinforce the importance of performing self-consistent solar model calcula¬ 
tions to account for the effect of exotic particles. Also, they show that assuming a universality 
in the constraints that are imposed, e.g. a given Li/ L 0 value, can lead to biased bounds on 
the properties of particles. 
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Figure 4: Relative changes of the neutrino fluxes with respect to SSM prediction as function of the 
luminosity contribution of axions (red line) and hidden photons (blue line). 


4.2 Solar sound speed profile 

In Figure 5 we represent the sound speed profiles (<5c s = Cg '°) )S Cs - th ) Q f models including ax- 

y Cs,obs J 

ions (left panel) and hidden photons (right panel). For each case, five different values of <710 
and x rn are considered. SSMs correspond to the null value of each parameter. The solid lines 
are obtained for solar models implementing AGSS09 surface composition while the dotted 
lines correspond to GS98 surface admixture. The shaded zones show the theoretical uncer¬ 
tainties (blue and red colors; not including composition uncertainties) and the uncertainties 
coming from the inversion technique (grey color). 

We see that axions produce effects on the sound speed in the inner region of the Sun. 
This is in line with the results obtained for the convective radius Rcz> which showed little 
variation in the presence of axion energy losses. This implies that the constraints on <710 
would be obtained from the measurements of the sound speed at t/Rq < 0.35. On the 
other hand, hidden photons produce noticeable effect also in more external regions, meaning 
that the entire sound speed profile could potentially contribute to constrain the product x m - 
However, this is not so simple because in our statistical method composition is free to vary. 
Variations in the sound speed profile will be compensated, at least partially, by changes in 
the solar composition. 

As it is seen from Figure 5, solar models implementing AGSS09 surface admixture 
provide a poor description of the sound speed profile inferred from helioseismic data. In 
particular, the sound speed prediction deviates at the bottom of the convective envelope by 
about ~ 1% with the helioseismic values. In this region, the modifications introduced either 
by axion or hidden photon energy losses are generally small and cannot explain the observed 
discrepancy. On the contrary, for t/Rq < 0.35, where axion or hidden photons effects are 
most relevant, the sound speed profiles of models with fixed reference surface composition 
deviate in the opposite direction as would be required to resolve the abundance problem. 

Although the composition is kept fixed in models shown in Figure 5, results suggest that 
the information encoded in the solar sound speed profile should be able to help constraining 
the axion and hidden photon properties. 
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Figure 5: Sound speed profile of models with axions (left panel) and hidden photons (right panels) 
for different values of the axion-photon coupling constant g±o and of the product x'm (expressed in 
10~ 12 eV) for hidden photons. Models are calibrated to the reference solar compositions GS98 (dashed 
lines) and AGSS09 (solid lines). Red and blue shaded zone corresponds to the model errors and the 
grey one to errors in the helioseismic inversion. 
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Figure 6: Comparison of sound speed profiles of SSMs. Red and blue lines are SSMs with AGSS09 
and GS98 reference compositions and all input SSM parameters fixed to their central values. The 
black line shows results for the best SSM [37] resulting from finding the SSM with free composition 
that satisfies Eq. 3.3. 


4.3 Best standard solar model 

Our statistical approach derives from that in [37], where it was applied to SSMs. By letting 
the solar composition free and and by adjusting the input parameters in SSMs (nuclear cross 
sections, microsocopic diffusion rate, etc.) thorugh the pulls £/ (that contributes to the total 
X 2 through the penalty factor it is possible to find the SSM that best reproduces 

helioseismic and solar neutrino data. Details on the calculations and results are given in [37]. 
Here, we just present a quick overview. In Figure 6 we compare the sound speed profiles of 
the two SSMs with GS98 and AGSS09 solar compositions and the best SSM (best fit SSM) 
resulting from [37]. 
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This figure illustrates the improvements of the best fit SSM over both the AGSS09 
and even the GS98 SSM that describes heliseismic properties quite well. In particular, the 
peak in the sound speed profile right below the convective zone at rj Rq = 0.713 is now 
smaller than 0.15% and the behavior of the sound speed in intermediate regions between 
r/R .0 = 0.5 and 0.2 is also greatly improved. The best fit SSM is achieved with composition 
parameters 5z vo \ = 0.45 and Sz re f = 0.19 that correspond to the logarithmic astronomical 
abundances e o = 8.85 and £pe = 7.52 for oxygen and iron respectively (the most relevant 
elements among volatiles and refractories, respectively), very close to GS98 values. The 
systematic pulls generally indicate variations at a 1 -a level or less from the central values for 
the input parameters of solar models. The most relevant systematic pulls of the best fit SSM 
are: £ opa = 1.07, £ d iff = -0.50, £ S34 = -1.2, £si 7 = -1-30. For this model, x[ 2 nin = 38.5 for 34 
observational constraints (see Sect. 2.4) indicating a good agreement between the model and 
the data. For comparison, the y 2 of the SSM is 179.6 when the composition is kept fixed at 
the AGSS09 values and 62.4 when fixed to GS98. 

The best fit model is almost completely independent on the reference solar composition 
that is initially used. This is partly due to the fact that the coefficients Bqj (eq. 2.7) are 
model independent to a very good approximation. Also, because variations between AGSS09 
and GS98 compositions are very similar for all elements grouped as volatiles (~0.14dex) and 
then for those grouped as refractories (~0.05dex). Therefore, the two-parameter analysis 
performed in [37], and upon which we base the analysis in this paper, can be considered as a 
family of solar compositions that smoothly connects the two GS98 and AGSS09 SSMs. The 
reference solar composition, then, becomes irrelevant for our analysis because the problems 
associated to the solar abundance problem can be efficiently compensated by the varying the 
composition and the systematic pulls. 

The best fit SSM represents our best effort to reproduce helioseismic and solar neutrino 
data. Solar models constructed to reproduce seismic data have been obtained previously by 
other authors. A couple of examples are the so-called solar seismic model [80] and the solar 
model by [52]. However, those models are not fully consistent with evolutionary history of 
the Sun because they are constructed as snapshots that reproduce present day structure. 
Our approach reproduces the present day structure consistently with evolutionary history as 
well. 

As a final comment, it should be added that the description of the solar properties 
that the best fit SSM offers should be understood as a phenomenological description of the 
solar interior structure. Alternatively, one could let the opacity profile free and constrain the 
composition (e.g. to AGSS09 values) and recover the same level of agreement between solar 
data and models [37, 63, 79]. 

4.4 Global results 

In this section we present the results of our global analysis. The bounds on axions and 
hidden photons are obtained by marginalizing with respect to surface composition, i.e. for 
each assumed value of for gio and x m we rescale the surface abundances of volatile and 
refractory elements by the factors (1 + 5z v Q i) and (1 + fc me t) in order to achieve the best 
possible agreement with observational data. The null values of the chemical composition 
parameters 5z vo \ and 6z me t correspond to the GS98 or AGSS09 photospheric composition 
(depending on the composition used to calculate the reference solar model). As expressed 
before, results have a very minimal dependence on the reference solar composition used. For 
simplicity, we show here the results obtained by using the AGSS09 as reference composition 
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Figure 7: Left panel : values of Na and Ay 2 for models with axions. Solid line: using all observables 
4>( 7 Be), 4>( 8 B), Y s , Rcz and 30 points of the sound speed profile. Dashed line: using the sound 
speed. Dotted-dashed line: using the neutrinos and convective envelope properties. Right panel: best 
fit composition as a function of <?io presented as the logarithmic astronomical abundances (red solid 
lines). Black lines represent the GS98 (dashed) and AGSS09 (solid) values for eo and tFe- 


(i.e. as pivot point for expansion in 5z re f and Sz vo \). Identical results are obtained if GS98 
composition is instead used. 

Axions - There is no model including axions (g \o ^ 0) that improves the overall fit to the 
data with respect to the our best fit SSM (Sect. 4.3). The variation of x 2 and its equivalent 
Na = v 7 Ax 2 as a function of gio are shown with solid line in left panel of Fig. 7. The right 
panel shows the values of the logarithmic abundances eo and eF e that provides the best fit 
to the data as a function of < 710 - These quantities are almost independent on g\o indicating 
that there are no degeneracies between composition and axion effects. Moreover, the best fit 
composition is very close to GS98 values. This should not be intended as a proof that the 
AGSS09 admixture is wrong but more as an evidence that the thermal stratification of SSM 
implementing GS98 composition is a good approximation of that of the real Sun. 

By setting a limit at Ay 2 = 9, we derive the upper bound 510 < 4.1 at a 3 -er CL. This 
is almost a factor of 2 lower than previous solar limits (Sect. 5). 

One of the advantages on using equation 3.2 to calculate y 2 is the possibility to under¬ 
stand which are the contributions to the total value of y 2 . As it has been discussed in [78] 
and [37], the method allows us to separate y 2 bs and y 2 ys . If we compare those values for the 
SSM and for the model with g\o = 4, we find that for the SSM y 2 bs = 33.6 and y 2 ys = 4.9 
while for the model with g\o = 4 we obtain y 2 bs = 35.3 and y 2 ys = 11.7. The values for 
y 2 bs are similar for both cases. This is the result of letting the composition to adjust to 
provide a good fit the solar data. Also, the changes that axions induce in solar models (see 
Sects. 4.1 and 4.2) are partially compensated by the systematic pulls. But changes to the 
input parameters come at the expense of increasing y 2 ys , that thus provides the dominant 
contribution to Ay 2 . The dominant pulls that increase the value of the x 2 with increasing 
gio are the solar luminosity (decrease), Sn (increase) and, to a lesser extent, S 17 (decrease). 
The first two are mostly related to changes in the sound speed in the solar core. A lower 
solar luminosity and a larger Sn both contribute to decrease the theoretical sound speed in 
the solar core [81], compensating the effect of axions (Fig. 5). Changes in S 17 occur in order 
to limit the increase in <L( 8 B) shown in Fig. 2. 

It is also instructive to discuss how much each piece of experimental information con¬ 
tributes to the bound on g±o . To this aim, we show in the left panel of Fig. 7 with dotted lines 
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the Ay 2 functions obtained by considering separately the sound speed profile on one hand, 
and the neutrino fluxes <h( s B) and 4>( 7 Be) and convective envelope properties Yg and Rcz 
on the other. We refer to this last dataset with the label i'+conv in the following. The most 
restrictive limit comes from the sound speed profile that, alone, sets the limit g±o < 4.6 at 3a 
CL. We note, however, that also the zz+conv dataset provides a restrictive bound g\o < 5.5 
at 3a CL. The latter value is more restrictive than found [36] using only the <L( 8 B) even if 
composition is free in our fit. This result that may seem surprising can be explained from the 
results presented in [37]. It was indeed shown in that paper that the two observable quantities 
Yg and Rcz permit to determine the surface composition of the Sun (in the two parameter 
analysis in terms of Sz vo \ and <5~ m et) with good accuracy. The experimental information on 
the neutrino fluxes, which are strongly dependent on < 710 , can be then effectively translated 
into a bound on axion energy losses. 
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Figure 8 : The values of Na and Ax 2 for models with axions using $( 7 Be), 4>( 8 B), Y s , Rcz and 
20 points of the sound speed profile as input observables. 

Finally, since the sound speed profile provide the most restrictive constraint, we check 
that the bump of c s (r) observed in the tachocline, i.e. the region just below the convective 
boundary, does not affect the final results in a critical way. This bump is present even in 
models with optimized composition (see e.g. the sound speed profile of the best fit SSM in 
Fig. 6) and is due to the inadequacies of SSMs in modeling dynamic effects in that region that 
affect the composition gradient and the transition between adiabatic and radiative gradient, 
(e.g. [82-84]). We repeated the global analysis excluding the sound speed data profile in the 
region r/R q > 0.6, thus reducing the sound speed data set to 20 data points. The exclusion 
of the sound speed determinations in the external radiative region, as it expected, improves 
the quality of the fit being Xmin = 12.4 for 24 observational constraints. However, since the 
limits are derived from the Ay ; 2 distribution, the bound for <710 does not significantly change. 
In fact, by excluding this regions, as can be seen in Figure 8 , g±o < 3.9 at 3a CL which is an 
even more stringent limit than the case where the full sound speed profile (in combination 
with the other experimental informations) is used. The fact that axions mainly affect the 
sound speed in the inner region of the Sun (see Figure 5) explains why the results depend 
only weakly on the inclusion of the bump. And, using the full sound speed the bound of g±Q 
is more conservative. 

Considering results discussed in this section, our recommended upper limit is g ai < 
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4.1 • 10~ 10 GeV^ 1 at 3-a CL for the axion-photon coupling constant. In Fig. 9 we summarize 
the relevant astrophysical constraints for hadronic KSVZ axions and axion-like particles, 
including our newly derived limit, together with the prospects of the future IAXO. 



Figure 9: Constraints on axion-like particles with a two-photon coupling as a function of the mass. 
The hadronic axion would be a point in the KSVZ line. The global fit presented in this paper 
provides a constraint (Sun (global)) that improves over the solar constraint obtained only with the 
solar neutrino argument (Sun v) and previous seismology constraints (Sun c s ) by a significant factor. 
Still, more stringent bounds come from the lifetime of Horizontal Branch stars in globular clusters. 
For masses below ~eV the axion heliosocopes CAST and SUMICO provide very competitive limits 
and the future International Axion Observatory (IAXO) has the potential of improving over them all. 
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Figure 10: Same than figure but for models with hidden photons. 
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Hidden photons - In Figure 10, we show the results of our analysis for models including 
hidden photons. The left panel shows, the Ay 2 distribution and the corresponding Na values 
as a function of the product \m, as obtained by using all observational constraints. The right 
panel shows the values of so and epe that minimize the y 2 for each assumed values of x m - I 11 
this case, as for the axion models, the best fit corresponds to the SSM with x m = 0 and the 
volatile and refractory abundances increased by 5z vo \ = 0.45 and Sz Te f = 0.19 with respect to 
AGSS09 composition that translate to £o = 8.85 and £Fe = 7.52. The bound at the 3<r CL 
is given by x m < 1-8 • 10~ 12 eV when the complete sound speed profile is used. If the region 
r > 0.6 Rq is excluded, then the limit is only marginally different, x m < 1-7 • 10~ 12 eV. Here 
again, using the full sound speed profile gives a more conservative limit than excluding the 
sound speed bump at the base of the convective envelope. 

tThe different observational data contributes to the final result as in the axion case. The 
sound speed is the most restrictive observable, giving the constraint x m < 2.4- 10~ 12 eV at 3a 
CL. The neutrinos fluxes and convective envelope properties also give the relevant constraint 
Xm < 2.8 • 10~ 12 eV at 3a CL. 

As a final result of our analysis, we quote x m A 1.8 • lCC 12 eV at 3-cr CL as an upper 
bound to product of the kinetic mixing parameter and mass of hidden photons. This limit 
does not depend on the assumed solar surface composition. In Fig. 11 we discuss our new 
constraint in the mass-mixing plane together with other concurrent limits. Notably, the 
constraint is the most stringent in the range from 3 • ICC 5 eV to 8 eV. 
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Figure 11: Constraints on hidden photons with mass m and kinetic mixing with photons %. The 
global fit presented in this paper, Sun-L (global), provides a constraint that improves over the so¬ 
lar constraint obtained only with the solar neutrino argument, Sun-L (^), and the direct detection 
constraint from XENON10 data [85]. 
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5 Discussion 


Here we have revised the solar limits on the axion-photon coupling constant and the product 
of the kinetic mixing and mass of the hidden photons by combining different solar constraints. 
We have used a statistical approach that accounts for both experimental and theoretical errors 
and lets the composition free in order that the solar abundance problem do not bias our final 
result. 

Previous works on the subject have often relied on setting an upper limit to the total 
energy loss carried away from the Sun expressed as a fraction of the solar photon luminosity 
Lq. It is then useful to present our results as a function of Lj/L 0 . Also, L*/L 0 sets a natural 
scale for comparing the effects of axion and hidden photon losses to each other. 
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Figure 12: A% 2 and Na as function of the luminosity contribution of axions (red line) and hidden 
photons (blue line). 

In Figure 12, global results from the previous section are presented again (only for the 
AGSS09 sets of models) but now as a function of Lj/L 0 . We observe that the 3-cr constraint 
for hidden photons corresponds to a contribution on the solar luminosity of only about 2% 
while for the axions this value is around 3%. For comparison, in [36] the upper limit to axions 
is set by demanding L a < O.1L 0 . We find that a much more restrictive limit can actually be 
imposed. 

It is interesting to note that a smaller fraction of L 0 is required for the hidden photons 
than for axions to set a given confidence level. We can explain this by considering how the 
observables change in both cases. In Fig. 5 we have seen that the changes on the solar sound 
speed profile are localized towards the center for the axion case while for the hidden photons 
the whole profile is affected. Also, we find a steeper relation between the energy lost through 
hidden photons and the change in <f>( 8 B) than in the case of axions. All this contributes to 
placing somewhat stronger constraints for hidden photons than for axions in terms of the 
fraction of energy that is lost through these channels. 

Solar limits on the coupling of hadronic axions to photons are not the most restrictive 
ones. However, this is a well-studied subject and it allows us to gauge the performance of our 
global statistical approach by comparing our findings to results from previous investigations. 

The first work that placed a limit on the axion-photon coupling constant using he¬ 
lioseismology is [35], where a set of solar models was calculated including effects of axions 
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self-consistently. The bound gio < 10 was found by restricting the axion luminosity at 
L a < O.2L 0 . For this limit, the deviations the sound speed profile are 0.8% at R = 0.1 Rq. 

More recently, [36] found g io < 7 by using the experimental <L( 8 B) determination and, 
based upon its agreement with SSMs values, assuming that the actual solar flux cannot 
exceed the SSM prediction by more than 50% (roughly equivalent to a limit of 3<r model 
uncertainty). As mentioned before, this is comparable to setting the limit L a < O.1L 0 . 

A different approach has been followed in [52], where a seismic model is used to de¬ 
termine the effect of the axions. Seismic models are static models of the Sun which, by 
construction, reproduce the sound speed profile. They do not account for the evolutionary 
history which determines the chemical composition profile of the present Sun. Axions were 
included in the seismic models by adding the energy loss in the energy equation and modify¬ 
ing the central temperature and the abundances with respect the standard case so that the 
sound speed and solar luminosity are recovered 2 . The constraint on g\Q is derived from the 
experimental value of t f>( 8 B) and it is set to g±o < 2.5 at a ltr CL. To compare, from Figure 
10 we find a limit at a la CL around <710 < 1.5. 

In comparison to these works, our result 510 < 4.1 to a 3-cr CL is much more restricting. 
This is a direct consequence of our global approach of combining consistently helioseismic 
and solar neutrino constraints. As discussed in the previous section, the sound speed profile 
is the most restrictive observational constraint. If we restrict our analysis to using zz+conv , 
we obtain gio < 5.5 and a corresponding L a < 0.06 L 0 limit on the axion luminosity. 

Using the limits on the parameters we have found for axions and hidden photons, we 
present upper limits for the respective fluxes on Earth expected in direct detection exper¬ 
iments such as CAST or IAXO. Using Eq. 15 from [53] and our limit g\$ = 4.1 we obtain 
4> a ~ 6.0 • 10 12 cm 1 • s 1 . During the CAST data taking, the limiting flux has been found to 
be more restrictive for a wide range of axion masses. [53] find a limit for the axion-photon 
coupling constant of g\Q = 0.88 for m a < 0.2eV corresponding to 4> a ~ 2.9 • 10 11 cm 2 • s 1 
and in [ 86 ] they find an upper limit of g\$ = 2.17 for the mass range of 0.02 < m a < 0.39eV 
corresponding to a solar flux of 4> a ~ 1.8 • 10 12 cm ~ 2 • s . As already stated in the introduc¬ 
tion, our limit does not improve the best limits on g ai but we can confirm that, at the CAST 
limits, it is not expected that axions would have a measurable effect with helioseismology 
and solar neutrinos. 

Hidden photons have a younger history than axions in the literature, and previous 
bounds on the kinetic mixing parameter based on solar models are limited to [ 8 , 9]. In 
the first case, a very conservative limit was derived by assuming that L} ip < L 0 , leading to 
Xm < 1.4 • 10 _ 11 eV. In the second, the more restrictive upper bound x m < 4 • lCU 12 eV 
for masses smaller than rn < 0.3keV was derived from the condition Lh p < 0.1 L 0 and 
Equation 4.1 . We have improved this limit by including the sound speed profile in the analysis 
and a more consistent treatment of uncertainties in our global approach. As discussed above, 
the limit x rn < 1-8 • 10~ 12 eV is obtained from models for which < 0.02 L 0 , a much 
smaller fraction than employed in previous works. Using Eqs 4.10 and 4.11 from [9], the 
upper limit for the flux on the Earth we derive is 4 >hp ~ 3.2 • 10 14 cm ~ 2 • s ^ 1 corresponding 
to x rn < 1-8 • 10 _ 12 eV. 

We remark that our bounds on axion-photon coupling and hidden photons kinetic mixing 
from are not affected by the ongoing solar abundance problem, a crisis brought about by 
the latest generation of spectroscopic determinations of the solar photosphere composition. 

2 Note that, in this approach, it is not discussed whether the solar seismic model of the Sun with axions 
energy losses can be realized as the result of solar evolution. 
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Indeed, we systematically studied the role of the composition used in solar models for the 
constraints on g ai and x m by letting the composition free in the study and marginalizing it 
for each value of gio and x m - The results showed that the values for the composition are 
more or less constant with increasing g\o and x m and that the values are close to the GS98 
and the best fit SSM, implying that models with this composition are a good representation 
of the actual thermal stratification of the Sun. This is reflected in the values of y Q bs> that 
are almost invariant at the gio or x m 3-cj limits compared to the best fit SSM. 

As shown in the literature [37, 63, 79], an equally good best fit SSM could be achieved by 
letting the radiative opacity vary instead of the solar composition. Keeping this in mind one 
concludes then that our limits are, to a good approximation, independent of the composition. 
This also indicates that exotic energy losses, even assuming the broad radial distribution pro¬ 
duced by hidden photons, cannot be advocated as a possible explanation of the disagreement 
between solar models with AGSS09 composition and helioseismic data. 

There are different ways in which to move forward. Our global approach represents a 
qualitative step forward in combining different sources of solar data for using the Sun as a 
laboratory for particle physics. But further data is available which we have not considered. 
For example, frequencies of low-degree oscillations can be combined in an advantageous man¬ 
ner to gain further insight on the solar core (e.g. [87, 88]). Using these data can enhance the 
constraining power of helioseismology. The caveat is that correlations among different helio¬ 
seismic observables have to be taken into account. We have already acknowledged that even 
the radial profile of the sound speed presents correlations that have not yet been quantified 
in the literature. This is the next step we will pursue for improving our statistical approach. 

It is important to notice that current uncertainties are dominated by model uncertain¬ 
ties. Moreover, there is not a unique dominant error source. For example, the 11% error in 
the theoretical <h( 8 B) results from 8% error in Sn and 7% from radiative opacities. If one 
includes solar composition as a source of uncertainty, this contributes a further 9%. Thus, 
reducing modeling errors is a difficult task. A similar situation happens for the sound speed 
profile, where the solar composition and the radiative opacities play comparable roles. 

On the other hand, the &(pp) and &(pep) fluxes have the interesting property that 
their theoretical errors are small, 0.6% and 1.2% respectively and therefore experimental 
determinations would place constraints on non-standard energy losses almost completely 
independent on solar modeling uncertainties. Recently [89], Borexino has directly measured 
the solar &(pp) for the first time, with an observational error of ~ 10% and the solar &(pep) 
[90] with an observational error of ~ 20%. Due to these large experimental errors, we have not 
included these fluxes as observational inputs because they would not contribute to the derived 
limits. In a future generation of experiments, reaching a 1% error in such measurements could 
provide a cross comparison between the solar photon and neutrino luminosities. Considering 
the case of hidden photons, for which the solar constraint is the most restrictive one over a 
wide range of photon masses, such a measurement would translate into a hard limit x m < 
3 • 10~ 12 eV at 3 -a CL, provided the measured value agrees with SSM results, and smaller if 
it does not. This limit, albeit larger than our estimate, has the attractive advantage that is 
practically independent of solar modeling. 

6 Summary 

We have extended the statistical approach presented [37] that combines in a consistent man¬ 
ner helioseismic and solar neutrino data for using the Sun for particle physics studies. In 


- 22 


order to avoid that the solar abundance problem could bias the final result, we considered 
the surface abundances of the Sun as free parameters in our analysis. 

As a test case we have applied the method to the well studied case of hadronic axions 
and showed that previous solar bounds on the axion-photon coupling can be improved with 
our methodology. We have derived a strong upper limit of g ai = 4.1 • 10 -10 GeV at a 3a 
CL, almost a factor of two better than previous results. As a further application, we have 
considered energy losses by hidden photons. This case is particularly interesting because the 
Sun offers the most restrictive limits over a wide range of photon masses. We have derived a 
new upper limit for the product of the kinetic mixing parameter and mass x m = 1-8- 10 -12 eV 
for a hidden photon mass in the range m < 0.3 keV, more than a factor of 2 improvement 
over previous results [9] and better than the direct detection constraint of XENONIO [85]. 

By comparing in detail results from axions and hidden photons, we also conclude that 
relations between dark luminosities and observables such as <f>( 8 B), e.g. Equation 4.1, depend 
on the type of particle under consideration and should be employed consistently to avoid 
biasing results. Finally, we have checked that including hidden photons in solar models with 
AGSS09 composition degrades the agreement between models and helioseismology, and thus 
does not help in mitigating the solar abundance problem. 

The approach is of course general and will be extended to account for correlations among 
observables, a fact hitherto neglected in the literature, that will allow us to include additional 
helioseismic constraints. 
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